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Abstract 

We discuss the problem of mass, noting that meson masses decrease with increasing 
scale as the dynamically generated condensate of “soft glue” is melted (Brown/Rho 
scaling). We then extend the Bielefeld LGS (Lattice Gauge Simulation) color singlet 
interaction computed for heavy quarks in a model-dependent way by including the 
Ampere law velocity-velocity interaction. Parameterizing the resulting interaction 
in terms of effective strength of the potential and including screening, we hnd that 
the masses of vr, cr, p and Ai excitations, 32 degrees of freedom in all, go to zero (in 
the chiral limit) as T ^ Tc essentially independently of the input quark (thermal) 
masses in the range of 1 — 2 GeV, calculated also in Bielefeld. We discuss other LGS 
which show qq bound states, which we interpret as our chirally restored mesons, for 
T > Tc. 

Key words: Chiral symmetries, Relativistic heavy ion collisions. Lattice QCD 
calculations 
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1 Introduction 


The problem of mass is one of the most fundamental in physics. We know 
that mesons have masses, all but the pion (in the chiral limit) at low densities 
and temperatures; i.e., at low scales. We now have experimental evidence 
that meson masses decrease by ~ 20% as the density increases to nuclear 
matter density. This incipient decrease has been seen in the STAR (STAR 
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Fig. 1. The decrease in vector meson mass can be simply understood as arising 
from the decrease in scalar density with increasing density. 

Collaboration) data for the p-meson, at a low density ~ O.lSno, where no is 
nuclear matter density [1]. 

Data on in-medium lo photoproduction in iV6 measured by the CBELSA/TAPS 
collaboration shows this mass to decrease roughly consistently with about 15% 
for nuclear matter density [2]. Such decrease in mass predicted by Brown and 
Rho [3] was made obvious in the work referred to above as coming from tad¬ 
poles connecting the scalar density to the relevant vector meson, as shown in 
Fig. 1, basically an extension of Walecka mean held theory. 

This is well and hue for hnite density, although the experiments take us up to 
only ~ no. However, Harada and Yamawaki [4] carrying out a renormalization 
group calculation in their vector manifestation show that the vector mass goes 
to zero at a hxed point as the temperature goes up to Tc from below. 

One can understand the behaviour of meson masses with temperature in terms 
of the two scales for condensed glue. 

The hard glue, or “epoxy”, is known to break scale invariance and provide the 
parameter Xqcd- R is an explicit breaking by the hard glue condensate which 
exists up to Tc^quenched — 250 MeV (See Fig. 2). 

In addition to the hard glue, there is a “soft glue”, condensed below T^, which 
connects the quarks in mesons and provides their masses. The soft glue pro¬ 
vides a dynamical breaking of scale invariance, the latter being restored at Tc 
[3]. 

The temperature dependence of the soft glue has been calculated in LGS by 
David Miller [5] and we show it in Fig. 2. The connection of the soft glue 
with masses is most easily understood in NJL (Nambu Jona-Lasinio) [7]; it is 
the difference in energy obtained by giving the quarks masses in the negative 
energy sea. 


B.F.(soft glue) 
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where m* is the medium dependent quark mass. As the soft glue is “melted” 
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Fig. 2. Gluon condensates taken from Miller [5] (Fig. 2 of BGLR[6]). The lines show 
the trace anomaly for SU(3) denoted by the open circles for light quark masses and 
the heavier ones by filled circles. The Tc marked in the figure is that for quenched 
QCD, whereas Miller deals with unquenched QCD with Tc = 175 MeV. 


with increasing temperature, the quark mass m* —0 at chiral restoration, 
Tc — 175 MeV. (The last point on the right in Fig. 2 is hard glue, or epoxy.) 
The soft gluon condensate is a dynamical breaking of scale invariance. Note 
that there is essentially no change in the epoxy from the next to last point. 
The soft glue holds the constituent quarks together in mesons. 


The epoxy is melted only as T goes up to Tquenched ~ 270 MeV. In Shuryak’s 
model [8] it consists of instanton molecules which do not break chiral symme¬ 
try. 


Now the evaluation of the tadpole process of Fig. 1 as T goes up to Tc from 
below would involve unquenched LGS below Tc, which have not yet been 
carried out. However, the Bielefeld group has carried out LGS in full QGD 
above Tc- We report here on those calculations in pure gauge, where the group 
has sent us their results electronically. Their full QGD SU{2) x SU{2) results 
[9,10] are essentially the same as in pure glue except for rescaling with the 
appropriate Tc’s. 
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2 The qq Bound States at T 


In quenched calculations Petreczky et al. [11] find quark and gluon (thermal) 
masses^ of ~ 1.6 GeV and 1.4 GeV, respectively for T = l.ST^ where we take 
T);,quenched = 270 MeV. Thus, quarks and gluons are not the thermodynamic 
variables of QGD in this energy range; their Boltzmann factors would be 
negligible. 

On the other hand, the SU{2) x SU{2) LGS transition is second order and, 
as a result, the pion mass must remain zero, in the chiral limit, as T moves 
upwards through Tc. Furthermore, the pion and sigma mesons are equivalent 
excitations as chiral symmetry and helicity are restored at T^, so, in the chiral 
limit, these should be massless. We have no such guiding symmetry above 
for the (equivalent) vector and axial vector excitations, except that the vector 
mass goes to zero as T approaches Tc from below [4]. 

We have shown [12] that the color singlet interaction mixes the tt and p just 
above Tc so strongly that one cannot talk about linear excitations until the 
temperature is well above Tc, the width of the excitations being as large as 
the p mass, which was taken to be ~ 280 MeV as arising completely from the 
explicit chiral symmetry breaking. In any case any additional mass assigned 
to this very nonlinear excitation is kept small by the minimization of the 
Helmholtz free energy [12]. The heavy quark interaction in LGS accomplishes 
only a small part of the vr, p mixing. 

We have 32 degrees of freedom which all go nearly massless at Tc in the chirally 
restored mesons tt, a, p, Ai. This is sufficient for LGS to give the correct entropy 
at Tc- 


3 Lattice Gauge Calculations above Tc 


The Bielefeld group [13,14] have carried out lattice gauge calculations to obtain 
the heavy quark free energy for the region of temperatures above Tc. They have 
given us their results for pure gauge calculations electronically. These agree 
with their SU{2) x SU{2) full QGD calculations [10,15] in almost all respects, 
once a rescaling by the relevant temperatures is made. We discuss in this 
section their results for the color singlet (Goulomb) potentials. 

The hnite temperature analog of the static potential is not known and it 
is not clear whether at all it can be properly dehned. Therefore we use the 

^ These may be ~ 30% smaller in full QCD (in which they have not yet been 
calculated) if rescaled with the appropriate Tc. 
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internal energy as potential in the Klein-Gordon equation. The color singlet 
internal energy can be derived from the free energy obtained in lattice gauge 
calculations: 


^F^{r T) 

l/i(r,T)=Fi(r,T)-T^^ (2) 

Above Tc, at large distance, the dominant scale of the effective strength of 
QCD potential is the temperature, i.e. g ~ g{T) for rT 3> 1. However, at 
short distances, the effective strength is dominated by the hard processes and 
the dominant scale is the length, i.e. g ~ g{r). In order to incorporate these 
properties into the effective strength of the potential [16], motivated by the 
Debye screened perturbative results, we £t the color singlet free energy Fi{r, T) 
as follows [13], 


Sttt 

^'^'■’^^^log[(l/r + HT)/A] 

where A, B, and D are fitting parameters which depend on temperature. In 
principle A also depends on temperature, but for consistency, A = 232 MeV 
is used for all temperatures. ^ In the above expression (Eq. (3)) we subtracted 
the value of the free energy at r — oo, namely Fi{r = cx), T). This is because 
we are only interested in the binding energy. In Kaczmarek et ah [13] the 
authors note that “at sufficiently short distances the free energy agrees well 
with the zero temperature heavy quark potential and thus also leads to a 
temperature independent running coupling. The range of this short distance 
regime is temperature dependent and reduces from r ~ 0.5 fm at T ~ Tc to 
r ~ 0.03 fm at T ~ 12Tc.” Since the rms radius of our chirally restored mesons 
is 0.2 fm at Tc, it is clear that for the regime of T in the neighborhood of Tc, 
the running coupling is the same as the much studied zero temperature one. 

As in Brown et ah [17], in order to enforce the asymptotic freedom at the 
origin, we used the molecule radius as 


R ~ 


h 

2mq 


( 4 ) 


This matches quite well the short-distance lattice regularization. For the out¬ 
side region, r > R, the effective potential can be determined using Eqs. (2) and 
(3) by fitting the results of lattice gauge simulation summarized in Fig. 3. In 

^ Since we fitted the Bielefeld LGS data, final results are not sensitive to the choice 
of different parameterizations. 
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Fig. 3. The potential for various temperature(mq = lAGeV). Filled symbols denote 
the results of lattice gauge simulation [18] and solid lines are fitted curves. The thick 
solid line is for the color-Coulomb interaction with constant coupling as = 0.5, as 
used in [17]. 

order to extend the effective potential inside the molecnle radins, the effective 
strength of the potential af^ can be defined by 


(t T) 

V(r, T) = - ^ ^ ^ for r > R. (5) 

r 

Inside the molecnlar radins, following Brown et al. [17], one can introdnce the 
effective potential as 


= forr<fi, (6) 

which is the same as the Conlomb potential for a nniform charge distribntion. 

Lattice calcnlations from the Bielefeld gronp [13] are the interactions between 
heavy quarks, so we need to add (model dependent) interactions so as to use 
them for the light quarks making up our mesons. Brown [19] showed that in 
a stationary state the E.M. interaction Hamiltonian between fermions is 


Hint — —(1 — ai ■ 52), 
r 


( 7 ) 
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where di ^2 are velocity operators. In this paper we restrict our consideration 
to T ~ Tc- It will turn out that the mesons are essentially massless in this 
region so ai ■ a 2 = ±1- Note that a is essentially the helicity. Hence we expect 
for the light quark potential to be 



Tc) 


2VCvy for tti ■ 02 = -1 
0 for oi ■ 02 = +1 


since helicity is either ±1 near T^. 


( 8 ) 


As we shall see from the results of our calculations the factor of 2 introduced 
by the magnetic interactions is precisely the factor needed to make the (zero) 
masses of tt and a continuous across Tc, as dictated by chiral symmetry (since 
the pion mass is protected and the a is degenerate with the pion at Tc). 

Note that the situation above Tc with quarks having thermal masses is unfamil¬ 
iar to many research workers, but it is completely analogous to that in heavy 
atoms where the thermal mass is simply a (repulsive) Coulomb potential, 
which transforms like the 4th component of a 4-vector. Indeed, the magnetic 
interaction in large Z atoms has been studied by Brenner and Brown [20]. Of 
course both Coulomb and magnetic interactions are repulsive here. Following 
their Eq. (23) it is pointed out that the magnetic interaction is 60% of the 
Coulomb for Z = 100. In atomic physics we do not have chiral invariance to 
constrain our result, nor can we go beyond Za = 1 with Dirac wave func¬ 
tions from pointlike charges, although it can go further with our modification 
Eq. (4), but the physics is the same since we deal with Coulomb interactions 
and the results are in the same ballpark. Furthermore, the Brenner and Brown 
results £t the measured iF-absorption edges of heavy atoms. 

Using the effective potential enhanced by velocity-velocity interaction, Eq. (8), 
we have estimated the color-Coulomb binding energies and the size of bound 
state in Table 1. In these estimates, we introduced the thermal quark mass[17,21] 
in the range of 1 — 2 GeV. This mass enters in the equation of motion as in 
Weldon [21] which is different from the massless Dirac equation. The final 
equation is close to the Klein-Gordon equation within 10% error, so we solved 
the Klein-Gordon equation instead of the full dispersion relation. 

As noted in BLRS [17] we solved the Klein-Gordon equation in the form given 
by Hund and Pilkuhn [22] 


(e 





= 0 


(9) 


where p is the momentum operator, and the reduced energy and mass e = 
— mf — m\)/2E, /x = 17111712 /E with mi = m 2 = mg. 
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Table 1 

Binding Energy due to color-Coulomb interaction AEc with velocity-velocity inter¬ 
action OL\-a .2 = —1, 4-point interaction energy AE 4 = G|'0(O)P, the radius of bound 
state (r^), reduced 4-point interaction energy A£' 4 ^r, and the mass of bound state 
Mb.s.- Energies and masses are in GeV, and the radius is in fm. 


T/Te 

AEc 

AE 4 


AE 43 

Mb.s. 



rrig = 

l.OGeV 



1.07 

1.55 

10.74 

0.12 

0.039 

0.42 

1.13 

1.33 

5.56 

0.16 

0.050 

0.62 

1.18 

1.12 

3.13 

0.21 

0.059 

0.82 

1.40 

0.70 

1.35 

0.29 

0.070 

1.23 

1.64 

0.56 

1.19 

0.31 

0.075 

1.37 

1.95 

0.34 

0.73 

0.37 

0.082 

1.57 

2.61 

0.16 

0.53 

0.47 

0.122 

1.72 



rrig = 

1.4GeF 



1.07 

2.03 

21.30 

0.10 

0.045 

0.72 

1.13 

1.69 

10.52 

0.14 

0.057 

1.05 

1.18 

1.37 

5.66 

0.18 

0.065 

1.36 

1.40 

0.85 

2.35 

0.24 

0.071 

1.88 

1.64 

0.69 

2.11 

0.25 

0.074 

2.04 

1.95 

0.45 

1.26 

0.30 

0.076 

2.28 

2.61 

0.26 

1.04 

0.35 

0.095 

2.45 



rrig = 

2.0GeF 



1.07 

2.68 

44.71 

0.08 

0.051 

1.27 

1.13 

2.14 

21.03 

0.11 

0.064 

1.80 

1.18 

1.68 

10.82 

0.15 

0.073 

2.25 

1.40 

1.02 

4.24 

0.20 

0.075 

2.91 

1.64 

0.85 

3.80 

0.21 

0.075 

3.08 

1.95 

0.57 

2.22 

0.25 

0.073 

3.36 

2.61 

0.38 

1.94 

0.27 

0.081 

3.54 


In addition to the color-Coulomb interaction, there exist 4-point interactions 
due to the instanton molecule interaction as indicated by the gluon condensate 
summarized in Fig. 2 [5]. In Brown et al. [6], we found that the 4-point coupling 




Fig. 4. Mass of bound states. The fitting curves show that the mass of bound states 
approaches to zero as T goes to T^. 

constant to be G = 3.83 GeV“^ above T^. In onr previons approach [17], we 
estimated that the 4-point interaction, 

AE, = Gmo)\^, ( 10 ) 

is mnch stronger than the color-Conlomb interaction. These properties are 
consistent with onr current calculation as summarized in Table 1. However, 
since the size of Coulomb bound state is now much smaller (see Table 1) than 
the instanton size, the volume effect has to be properly considered. 

In this paper, we propose a new approach to the instanton molecule interac¬ 
tion. The instanton molecule at Tc results from the close packing of instanton 
and antiinstanton in the time direction, each of diameter 2/3 fm, so the length 
of the instanton molecule is 

1 4 

— = -fm. (11) 

Since the instanton radius is 1/3 fm, the effective volume of the instanton 
molecule is 




-fm X TT 
3 



0.465 fm^. 


( 12 ) 


By considering the hnite size of the Coulomb bound state, the 4-point inter¬ 
action can be reduced as 
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P((») p((0) 



Fig. 5. Spectral functions of Asakawa et al.[24]. Left panel: for N-j- = 54 (T ~ 1.4Tc). 
Right panel: for Nr = 40 (T ~ 1.9Tc). 





(13) 


The reduced 4-point interactions are summarized in Table 1. Note that the 
binding is now dominated by the color Coulomb interaction for smaller tem¬ 
peratures due to the smaller size of the bound state. The hnal masses of bound 
states are summarized in Table 1 and Fig. 4 for different thermal quark masses. 
The hgure clearly shows that the mass of bound state approaches zero as T 
decreases to T^, independently of the thermal quark masses. The bound states 
in Table 1 include the effects of the velocity-velocity interaction. 


We have used Eq. (8) in calculating Table 1, so our numbers are only meaning¬ 
ful at and slightly above T^. Higher up we should switch over to the velocity- 
velocity interaction of Shuryak and Zahed [23]. We know that at a higher 
temperature T = T^em binding the velocities will go to zero as the molecules 
break up. This is the region of the “perfect liquid” [24]. Our Table 1 is still 
useful for T = lAT^ where LGS have found our chirally restored mesons using 
the maximum entropy method [24]. 

In Table 2, we summarized the results with only the heavy quark interactions. 
One can interpolate between cki • 02 = 1 and zero, using the Shuryak and 
Zahed [23] in order to determine the velocities. 


4 Lattice Evidence for the qq Bound States above 


In Fig. 5 we show the spectral functions of Asakawa et al [24]. Similar re¬ 
sults were obtained by Petreczky [25]. Only the peak lowest in energy at each 
temperature should be considered. The higher ones are lattice artefacts. We 
especially consider the 2 GeV peak at T = l.dT^. The lattice calculations are 
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Table 2 

Binding Energy due to color-Coulomb interaction AEc with only heavy quark in¬ 
teraction (without velocity-velocity interaction <3 • 02 )• Notations are the same as 
in Table 1. 


T/T, 

AEc 

AE 4 

yn 


Mb.s. 



ruq = 

l.OGeP 



1.07 

0.87 

1.53 

0.27 

0.065 

1.06 

1.13 

0.56 

0.77 

0.34 

0.067 

1.38 

1.18 

0.37 

0.43 

0.42 

0.069 

1.56 

1.40 

0.12 

0.18 

0.62 

0.090 

1.79 

1.64 

0.04 

0.12 

0.87 

0.116 

1.84 

1.95 

0.002 

0.02 

3.27 

0.023 

1.98 

2.61 

- 

- 

- 

- 

- 



ruq = 

lAGeV 



1.07 

1.03 

2.60 

0.23 

0.071 

1.70 

1.13 

0.66 

1.27 

0.29 

0.069 

2.07 

1.18 

0.45 

0.69 

0.36 

0.067 

2.28 

1.40 

0.18 

0.31 

0.49 

0.075 

2.55 

1.64 

0.09 

0.25 

0.58 

0.102 

2.61 

1.95 

0.03 

0.12 

0.90 

0.116 

2.66 

2.61 

- 

- 

- 

- 

- 



mq = 

2.0GeV 



1.07 

1.22 

4.73 

0.20 

0.077 

2.70 

1.13 

0.78 

2.20 

0.25 

0.072 

3.15 

1.18 

0.54 

1.17 

0.30 

0.068 

3.39 

1.40 

0.24 

0.54 

0.39 

0.067 

3.69 

1.64 

0.15 

0.48 

0.42 

0.078 

3.77 

1.95 

0.07 

0.28 

0.54 

0.094 

3.83 

2.61 

0.01 

0.12 

1.17 

0.118 

3.87 


quenched, so here is ~ 270 MeV. Possibly all scales would come down by 
the factor (Tc,unquenched777,quenched) h fuH QCD could be Carried out. 

We should remember that these vibrations were calculated in the heavy quark 
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approximation. 


First of all, we note that there is a complete SU(4) symmetry. In BLRS[17] 
we found no dependence of masses on isospin and negligible dependence upon 
spin. Thus the degeneracy of the S, PS, V, and AV degrees of freedom should 
not be surprising. 

As noted in BLR [12], for temperatures not far above Tc, the color singlet 
interaction mixes the p and tt excitations strongly, so the linearity that is 
present in the T = 1.4Tc results would not be expected to be present at lower 
temperatures, even in the heavy quark approximation in which the effective 
color singlet coupling is only half that for light quarks in our approximation 
(and the square of the coupling, which comes into the nonlinearities, only 
one-fourth). 

In this sense, it is useful that the heavy quark interaction is small enough so 
that the vibrations are roughly linear at T = 1.4Tc when calculated with it. 

In any case, we see that the bound state masses increase rapidly as T moves 
upwards from Tc. However, Shuryak and Zahed [23] show that very large num¬ 
bers of other excitations, mostly colored, come in in such a way as to furnish 
the pressure which is lost by the massless excitations at Tc contracting out, to 
a large extent, as thermodynamic variables as their mass goes up with increase 
in temperature. We will not deal with the many further excitations above Tc, 
considered by Shuryak and Zahed [23], but simply note again that the increase 
in entropy in going up to Tc can be furnished by our 32 massless mesons there. 


5 Discussion 


In BLRS[17] we based our construction of the tt and a mesons on their chiral 
properties; i.e., that the zero mass (in the chiral limit) of the pion should go 
smoothly through T^. At Tc and above the a is degenerate with the pion. The 
symmetry in isospin followed from our Lagrangian in which r = (r, 1) is a 
4-vector, so that the p and uj were components of the same particle, also the 
a and tt. There was, however, no symmetry requirement on the mass of the 
p, as T went down to Tc from above. (In the vector manifestation of Harada 
and Yamawaki [4] m* must go to zero at a hxed point as T goes up to Tc from 
below.) 

Spin effects are, however, small. The magnetic moment of the quark or anti¬ 
quark is [17] 
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(14) 


fJ'q,q 4 “ 

Po 

where po ~ ^rriq^q and mq^q is the thermal mass. For our chirally restored 
mesons at Tc, we see that ~ 2, which with Eq. (23) of Brown et ah [17], 
would give nip ~ mq/6. 

Although the confluence of curves giving the mass of the bound state in Fig. 4 
looks very impressive at hrst sight, we caution that the data points (points 
from the LGS) that we have for the lower T/Tc values do not go to distances 
shorter than ~ 0.25 fm, whereas our rms radius at Tc is only ~ 0.5 of this. 
However, Kaczmarek et ah (denoted as KKZP below) [13] have shown that at 
sufficiently short distances, up to r = 0.5 fm at T ~ T^, the heavy quark free 
energy agrees well with the zero temperature heavy quark potential, which 
is essentially Coulombic. KKZP emphasize that in this distance range for Tc 
(and shorter ones for higher temperatures) processes in the QCD plasma phase 
are still dominated by properties of the QCD vacuum. The KKZP analysis 
suggests that “it is more appropriate to characterize the non-perturbative 
properties of the QCD plasma phase close to in terms of remnants of the 
conhnement part of the QCD force rather than a strong Coulombic force.” 

It’s hard to separate the pure Coulombic effect from non-perturbative effects 
in our approach. However, these properties are already included in our re¬ 
sults because we parameterized our potential to £t the lattice results. In our 
parameterization, these non-perturbative effects are absorbed in the effective 
strength of the potential, Eq. (5). We believe that the stronger couplings for 
lower temperatures towards indicate the stronger non-perturbative effect 
near Tc- For small distances, r ~ 1 fm, the strength of the potential 

g{r,T) is dominated by the radius, close to the zero temperature limit. 

We can obtain a semi quantitative picture of why the meson masses are small 
at Tc in the following way. KKZP show that the short distance color averaged 
free energy is dominated by the singlet contribution. This interaction is very 
large at large distances, as shown by the Polyakov loops in quenched QCD 
(Fig. 4 of KKZP [13]). Thus, overall we believe that we can talk about Coulomb 
interactions (realizing that the discussion will not be gauge invariant) in order 
to provide a simple picture. 

With the doubling of the effective Coulomb interaction by inclusion of the 
velocity-velocity term and with the Casimir operator 4/3, we would obtain a 
(large distance) 



ttefr(^^) 


47r 


16 

y 


from the Polyakov loops at Tc, corresponding to 


( 15 ) 
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QeS ~ 8. 


( 16 ) 


Of course this implies strong coupling. It seems reasonable that this large color 
singlet coupling is largely responsible for the large indeed, rescaling the 
quenched 1 GeV which we used at Tc, would give ~ 600 MeV for the 
unquenched one. 


Potential V 



Fig. 6. Schematic model of the color singlet potential and the meson binding. 

In order to show why all of the curves in Fig. 4 sweep to zero, independently 
of TUq, in Fig. 6, we make our picture in terms of Coulomb energies and po¬ 
tentials, two repulsive Coulomb energies 2mq being the meson energies at a 
large separation; i.e., we move our zero in energy down to —2mq. 

The Coulomb is attractive, but in our gauge theory, must go to zero as the 
separation r —0, because of asymptotic freedom. ^ 

Viewing the potentials as like that in the hydrogen atom, we would expect 
the bound level to come down about halfway in the Coulomb potential; i.e., 
the zero point energy would be about half of the potential depth. 

With the doubling in the potential from the velocity-velocity interaction, the 
molecular state comes down to —2mq. This argument explains why the states 
in Fig. 4 all go to zero, independently of TUq. Our calculations support this 
picture; the potentials are approximately rescaled hydrogenic ones. 

We see from our model how asymptotic freedom giving perturbative interac¬ 
tions at short distances, can be reconciled with strong large distance behavior. 
In fact, KKZP[13] hnd that for their purpose of dehning a running coupling 
at hnite temperature the singlet free energy is most appropriate, as it has at 
short (r ■ max(T, Tc) 1) as well as large (rT 1, T ^ Tc) distances and 
temperatures a simple asymptotic behavior which is dominated by one gluon 
exchange; i.e., 

^ As we noted earlier, asymptotic freedom is built into our schematic potentials 
[17]. 
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(17) 


Fi(r,T) 


r-max(T,T,)«l 

9^{T) -g{T)rT 
Sirr 


In our case the meson at Tc with rms radius y (r^) < 0.12 fm is well inside the 
~ 0.5 fm Debye screening length, so only the (upper) perturbative Coulombic 
expression in Fi is needed. 

Our essential step in going from the heavy-quark color singlet potential to the 
light-quark description is to add the velocity-velocity interaction. In BLRS[17] 
the instanton molecule interaction played an even larger part, but with our 
very much smaller Coulomb bound states here from the lattice parameters, 
these are essentially negligible. The velocity-velocity potential was derived for 
the interaction of K-electrons in heavy atoms, holding for stationary states 
[19]. This interaction was important in reproducing the K-absorption edges of 
heavy atoms [20], where v/c was not much less than unity. We have shown 
[12] that low-mass chirally restored mesons just above Tc will explain some 
of the data from RHIC, which is difficult to understand in the standard sce¬ 
nario, and have made predictions for other experiments. Chiefly, however, our 
introduction of the velocity-velocity dependent term allows us to separate the 
64 degrees of freedom into 32 that have zero energy at T^, in our model and 
in the chiral limit, at Tc, with the other 32 lying in mass ~ ^rug, and there¬ 
fore irrelevant for the thermodynamics. Our approach is the only one giving 
exactly the correct number (as measured in LGS) of degrees of freedom at Tc. 

Our results using the numbers from the Bielefeld LGS essentially conhrm 
the scenario of Brown et al. [26] that chiral restoration is reached by mesons 
going massless. This followed from Brown/Rho scaling [3]. The discussion of 
relevant degrees of freedom; i.e., relevant mesons was, however, clarified only 
much later in BLRS [17]. The 32 degrees of freedom going massless at Tc is 
precisely the number lattice calculations require in order to fit the entropy 
increase up to T^. The pressure is given by these degrees of freedom only for 
a narrow range of temperatures from Tc up to ~ 1.2Tc where contributions 
begin from an additional large number of binary bound states, both colorless 
and colored pairs gq, qq and gg [23]. 

Although chiral restoration was discussed many years ago in terms of chirally 
restored mesons, the importance of the color Coulomb interaction was only 
understood recently, when lattice calculations found charmonium to be bound 
up to ~ 2Tc. This statement does not take into account many earlier published 
studies by Ismail Zahed, who explained to us the importance of the color 
Coulomb interaction. We emphasize our indebtedness to his work and to his 
tutelage. 

In BLRS [17] we employed a schematic model for the heavy quark color 
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Coulomb interaction with as chosen to be ~ 0.5 in order to bind charmo- 
nium up to ~ 2Tc. Our Coulomb bound states were then rather large at Tc, 
with rms radius ~ 1/3 fm, roughly the size of the instanton. Given these large 
bound states, the instanton molecule interaction played an important role in 
bringing their energy to zero at Tc- 

In the neighborhood of Tc the lattice color Coulomb potential we use here 
is more than twice as large as the schematic one used in BLRS [17], with 
the results that the Coulomb bound states are so small that the instanton 
molecule interactions are negligible at Tc. 

Aside from reproducing correctly the number of degrees of freedom in the 
entropy, energy and pressure at Tc, our scenario should be able to give de¬ 
scriptions in detail of dynamical processes, such as the p°/7r“ ratio measured 
by STAR Collaboration [12], We also noted in that reference that many dilep¬ 
tons should be emitted from the region just above Tc where the p-mesons have 
low masses. In Brown et ah [12] it was incorrectly stated that the p-mesons 
from the region of temperatures just above Tc would give dileptons of low 
invariant masses, especially where the cocktail peaks up. In fact, because of 
the ~ 2.8T thermal energy of the p above Tc, the invariant dilepton masses, 
which contain this thermal energy, should come in the ~ 500 MeV region, sig- 
nihcantly increasing these from the lower end of the mixed phase or from the 
rhosobar, as seen in CERES. However, since the invariant dilepton mass comes 
almost completely from the thermal energy in the new RHIC component of 
dilepton, the p± distribution should be different from those from CERES, a 
much higher proportion of the new RHIC component having high p±. 

Our work addresses only the baryon number zero problem, since we use LGS 
results and we deal only with mesons. We believe that it is relevant for RHIC, 
in which after the baryons clear out early, the material is chiefly composed of 
mesons. 
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A Appendix: Trace Anomalies, etc. 


The movement towards zero mass of the 32 mesons as T approaches to Tc may 
appear somewhat miraculous, but the movement of the masses to zero as T 
goes up to Tc from below was anticipated by Brown et ah [26] although the 
strongly interacting form of matter above was not predicted; i.e., it was not 
known that the thermal quark masses would be so large, and that, therefore, 
interactions must be very strong to bring the meson masses to zero just above 
T 

Adami et ah [27] on the basis of Deng’s lattice calculations [28] found that 
about half of the glue melted below T^, and the other half above, as shown 
in our Fig. 2. Shuryak, in a number of papers, interpreted the glue condensed 
below Tc as arising chiefly from the random instantons, some of which were 
melted going up to Tc, the others transforming into instanton molecules which 
do not break chirality. We generally discuss the glue below Tc as soft glue, 
that above as epoxy, (hard glue). It can be seen in Fig. 2 that the glue above 
Tc at the last point on the right is at the same position as the point at T^, i.e., 
there has been no melting from Tc until further above Tc. 

Formally, it is more precise to talk about the dynamical breaking of scale 
invariance by the soft gluon condensate below Tc. This is the glue that is 
necessary in order to hold the quarks together to make hadrons, and as this 
glue melts the hadron masses go to zero as in Brown/Rho scaling [3]. 

The Yang-Mills equation are invariant under transformations of the full con¬ 
formal group. Here we discuss only scale transformation, the simplest of these. 
A simple model in terms of scalar helds was considered by Freund and Nambu 
[29]. A scale transformation (j){x) —> X(j){Xx) on a scalar held induces the 
transformation on the scale invariant part of the Lagrangian 


^ X Linv^Xx^. (A.l) 

Now the trace anomaly in QCD can be expressed in terms of the divergence 
of the dilaton current whose non zero value signals the breaking of scale 
invariance. Including the effect of quark masses, the trace anomaly is given by 


= ©C = ^ m,qq - (A.2) 

q 9 

where q and are the usual quark and gluon helds. The second main as¬ 
sumption, the one with which we are now concerned with, of Brown and Rho 
[3] is : “In order to be consistent with the scale property of QCD, effective La- 
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grangian must reproduce faithfully the trace anomaly A(l) in terms of effective 
fields. ” 

Important for us is the fact that there are two kinds of glue, soft and epoxy, and 
two kinds of scale breaking, dynamical (or spontaneous) by the soft glue, the 
glue that holds the quarks together to make hadrons, and the explicit breaking, 
which takes place already in quenched QCD, giving the gluon condensate the 
scale parameter Xqcd- 

Also important for Brown/Rho scaling is that the spontaneous breaking can¬ 
not occur without the explicit breaking [29]. Without the lattice, there is no 
preference to realize the dynamical breaking in one region rather than an¬ 
other, but once Tc for the unquenched system is established, it is natural that 
this gives the scale above which the soft glue has disappeared and only epoxy 
remains. 

Consequently, there are two scales, (1) Tfiunquenched): the scale at which 
the soft glue has been melted; more importantly, the temperature for chi¬ 
ral restoration, (2) Tfiquenched): the scale at which the hard glue has been 
melted. 

Brown/Rho scaling is understood most simply by nuclear physicists brought 
up in Walecka theory in terms of three-body interactions shown in Fig. 1, with 
scalar “tadpoles”. However, these are not easily evaluated in lattice gauge cal¬ 
culations, even at hnite temperature, because they would require unquenched 
calculations which are not yet accessible. Furthermore, even when they are 
carried out, the hrst ones will probably have quite large bare quark masses. 
On the other hand the calculation of the gluon condensate below Tc, shown 
in Fig. 2, depends only weakly on the bare quark mass, only through the 
ffqVnqqq term in the trace anomaly of Eq. (A.2), which is small compared 
with the term. Therefore, it is much simpler to study the way in 

which meson masses drop with temperature by lattice calculations of the soft 
glue. 

Of course, for the present paper concerned with T ^ Tc the soft glue is absent 
and that is why the mesons we calculate have no more scalar masses, and why 
their energies go to zero at Tc to make continuity with the massless tt and a 
just below Tc (the p and Ai following suit because of the negligible spin- and 
isospin-dependence). 

We have emphasized that our 32 chirally restored mesons give just the number 
of bosonic degrees of freedom obtained from the measurement of entropy at 
Tc in lattice calculations. Of course, colored excitations enter not far above 
Tc, contributing to the entropy which is nearly constant over a substantial 
increase in temperature. Whereas we refer the reader to Shuryak and Zahed 
[23] to understand this latter point, we wish here to emphasize the utility of 
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having a layer of colorless mesons just at Tc and slightly above. First of all, 
the interactions are very strong as emphasized in BLR [12], and essentially 
inelastic so that there is no doubt of equilibration, which is necessary in order 
to dehne a thermodynamic variable such as temperature. Secondly, how do 
the mesons know that they can get out of the hreball (be emitted) below Tc, 
whereas they are conhned above ? In our scenario they have to go through 
the colorless layer of mesons at Tc, having deposited their color above Tc as 
the temperature decreases with expansion of the hreball. In other words, our 
layer of colorless mesons at Tc is a “color puriher”. 

The way in which the mesons, massless just above Tc, acquire their on-shell 
mass as they emerge below Tc was only briehy discussed in BLR [12], mostly 
for the p mesons. Hadrons were not considered and the whole issue of how 
quantitatively the hadrons go back on shell will require further study. Analysis 
of the copious RHIC data should be useful here. 
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